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I.  INTRODUCTION 


The  performance  of  a  modern,  high  temperature  gas  turbine  engine  is 
compromised  significantly  by  associated  requirements  for  component  cooling 
This  cooling  is  normally  accomplished  with  air  bled  from  the  "cold", 
compressor  end  of  the  engine.  Consequently,  the  paths  provided  to  conduct 
this  cooling  air  require  careful  design  and  flow  analyses  for  effective 
utilization  of  bleed  air  resources.  Gas  path  leakage,  which  may  not  provide 
any  useful  function,  occurs  between  engine  components.  The  combined  cooling 
and  leakage  flows  must  be  determined  and  their  impact  on  engine  performance 
evaluated . 

The  determination  of  these  cooling  and  leakage  flows,  called  the  internal  flow 
analysis,  requires  the  mathematical  modeling  of  a  complex  network  of  conduits 


and  restrictions  located  inside  and  outside  the  main  gas  path  from  the  engine 
inlet  to  the  final  nozzle.  The  total  pressure  losses  through  these  conduits 
and  restrictions  must  be  characterized  so  that  the  flow  capacities  can  be 
calculated.  Two  parameters  are  used  somewhat  interchangeably  for  the  flow 
characteristics  or  the  total  pressure  loss  characteristics  of  the  constituent 
restrictions .  The  discharge  coefficient,  C^,  is  a  measure  of  the  flow 
passing  through  a  restriction  relative  to  the  calculated  ideal  flow  at  the 
actual  upstream  and  downstream  pressures.  The  total  pressure  loss 
coefficient,  k,  is  a  measure  of  the  energy  required  to  drive  the  actual 
flow  through  the  restriction.  For  example  at  the  same  operating 
conditions,  the  discharge  coefficient  and  the  total  pressure  loss 
coefficient  based  on  maximum  upstream  dynamic  pressure  are  related  as 


CD  =“ 


k  -  k  +1 
e 


for  a  restriction  in  the  incompressible  flow  regime. 


ll\lomenc lature  definition  is  found  on  page  122. 


raw 


The  purpose  of  this  program  is  to  formulate  flow  charac ter  1 st 1 c  models  for 
conduit  bends  and  turns,  branches,  sudden  area  changes,  and  orifices 
applicable  to  the  restrictions  in  the  internal  flow  systems  of  gas  turbine 
engines . 

Modeling  of  flow  through  the  interior  cooling  and  leakage  passages  of  gas 
turbine  engines  is  an  inexact  science.  These  passages  are  often  of 
unconventional  geometries  for  which  experimental  data  do  not  exist  and  for 
which  numerical  fluid  dynamic  analysis  is  unreliable  or  impractical.  As  the 
result  of  in  situ  performance  variations  due  to  local  conditions  of 
turbulence,  approaching  flow  profile  characteristics ,  proximity  of  downstream 
restrictions ,  heat  transfer,  and  engine-to-engine  configuration  and 
dimensional  variations,  careful  rig  tests  on  the  actual  engine  parts  will  not 
yield  precise  flow  characteri sties  for  the  internal  flow  system  model.  These 
uncontrol lable  consequences  of  gas  turbine  engine  design  distinguish  this  flow 
network  analysis  from  the  more  exact  solutions  for  conventional  piping  or 
ducting  systems.  The  accurate  modeling  of  internal  flow  systems  of  gas 
turbine  engines  now  relies  on  the  modification  of  "reference"  restriction 
characteristics  with  application-specific  empirical  factors  based  on  global 
experience  from  engine  testing.  These  limitations  do  not  preclude  the 
reasonable  preliminary  predictions  of  internal  flow  system  performance  for 
untested  engine  designs.  Fortunately,  the  internal  flow  system  is  typically 
comprised  of  many  restrictions  in  series  and  parallel  arrangements.  The 
composite  nature  of  such  flow  networks  generally  relegates  the  flow 
restriction  characteristics  to  secondary  importance  with  respect  to  the 
correct  evaluation  of  flow  areas. 

An  exhaustive  literature  search  indicated  that  the  k-factor  is  the  parameter 
of  choice  for  general  restriction  geometries.  However,  the  definition  of  the 
k-factor  does  not  enjoy  the  same  consensus.  The  flow  models  presented  in  this 
report  are  based  upon  total  pressure  loss  k-factors  which  are  referenced  to 
flow  conditions  calculated  at  the  minimum  cross-sectional  area  at  the  upstream 
end  of  the  restriction.  The  rationale  for  this  selection  is  discussed  in  the 
theoretical  analysis  section. 


The  important  objective!  of  k-factor  modeling  are: 

2 

1)  Prediction  of  realistic  trends  and  boundary  values  at  approx  1  mate  1 y 
correct  levels  of  k-factor  for  the  component  geometry  at  local  average 
flow  conditions . 

2)  Relatively  simple  (usable)  formulation  of  the  primary  geometric  and 
fluid  dynamic  parameters  into  a  correlation  representat i ve  of  a  broad 
range  of  configurations  and  flow  environments. 

Such  k-factor  models  allow  realistic  comparisons  for  the  evaluation  of  design 
changes  and  environment  modif ications  Prediction  of  absolute  performance 
levels  will  usually  require  an  experienced  adjustment  of  the  appropriate 
k-factors  to  match  experimental  results. 

Often  internal  flow  system  models  for  new  engine  designs  are  synthesized 
initially  with  k-factors  for  a  static  orifice,  k  =.  2.7  ,  and  for  an  lsentropic 
nozzle,  k  =  1  0  ,  in  conjunction  with  exact  calculations  of  the  controlling 


passage  areas  These  models  can  do  surprisingly  accurate  when  carefully 
formulated  by  an  experienced  flow  analyst  The  preliminary  internal  flow 
model  is  refined  using  component  k-factors  appropriate  to  the  engine  design 
details  Later,  when  engine  performance  testing  yields  measured  pressures  and 
temperatures  for  the  internal  cavities  and  passages,  the  k-factors  for  the 
controlling  restrictions  can  be  modified  to  simulate  the  in  situ  pressure 
changes 


Assessing  the  validity  and  accuracy  of  k-factor  models  for  even  basic 
restriction  geometries  is  difficult  or  impossible  without  extensive 
experimental  support.  Generalized  k-factor  models  are  nonexistent  for  the 


These  "boundary  values"  could  also  be  termed  limiting  or  extremum 
values  As  an  example,  the  flow  losses  for  bends  of  increasing  radius 
ratio,  r/4.,  or  decreasing  bend  angle,  0,  should  approach  the  wall 
friction  loss  for  a  straight  duct  of  equal  bend  length  and  the  same 
cross-sectional  geometry  as  a  lower  limit. 
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broader  range  of  rpstrution  qi'omr tr  i(>5  necessary  for  gas  turbine  interrw] 
f  lou  analysis  When  the  effects  of  unt  onuent  1  ona  1  installations  and  fluid 
flow  environments  are  considered,  the  ability  to  precisely  predict  restriction 
flow  c harac ter  1 s 1 1 c s  is  presently  unattainable  However,  the  purpose1  of  this 
study  is  to  develop  approximate  k-factor  models  for  generic  bends  and  turns, 
branches,  sudden  area  changes,  and  orifices  common  to  gas  turbine  secondary 
flowpaths  These  algorithms  will  produce  repre sentat l ve  trends  and  boundary 
values  for  tne  identifiable  variables  of  the  geometry  and  the  flow  processes 

Open  literature  contains  many  performance  models  for  the  basic  restrictions . 
Some  of  these  empirical  models  were  derived  from  poorly  controlled  or 
incompletely  formulated  experiments.  Consequently,  some  of  the  available 
k-factor  models  are  limited  to  unspecified  ranges  of  influential  geometric 
and/or  flow  parameters  A  few  k-factor  models  even  produce  physically 
inconsistent  performance  predictions  in  particular  operating  regimes  A  small 
sample  of  these  exceeding  limited  k-factor  predictions  are  derived  from 
oversimplified  analytical  models  of  the  flow  phenomena. 

Sometimes  recognized  expert  opinions  exist  about  the  reliability  of  certain 
restriction  models.  Beyond  this  the  only  viable  procedure  for  selecting  among 
the  potpourri  of  k-factor  models  must  rely  upon  comparisons  of  performance 
predictions  at  selected  conditions  for  several  of  the  more  comprehensive 
models  and  upon  evaluation  of  their  boundary  values  where  possible.  The 
development  of  the  k-factor  models,  or  perhaps  more  correctly  the  synthesis  of 
the  k-factor  models,  for  application  to  the  analysis  of  internal  flow  systems 
in  gas  turbine  engines  will  be  accomplished  with  such  a  procedure.  One  or 
more  sources  will  be  utilized  to  produce  a  consistent  algorithm  of  acceptable 
accuracy  that  will  predict  realistic  performance  trends  for  variations  in 
component  geometry  and  flow  conditions. 

The  k-factor  models  sought  in  this  study  will  be  formulated  with  influence 
coefficients  to  correct  the  "reference"  performance  predicted  for  a  basic 
geometry  and  flow  environment,  eg  , 

k  -  f(x,y,z,Re) 


A 


I'rri  :  .  -V*  i'v' 


-  -r-  -V 


for  the  effects  of  non-standard  geometry  or  unusual  flow  conditions. 


k  =  C.  C,  C  k* 

1  i  n 

Basic  flow  environments  are  generally  for  the  isolated  component  in  the 
incompressible  regime.  This  implies  fully  developed  entrance  flow  and  the 
effect  of  complete  downstream  pressure  recovery .  The  influence  coefficients 
for  variations  from  the  basic  geometry  or  standard  incompressible  flow 
characteri sties  will  be  developed  from  available  data  sources.  It  is  worth 
noting  that  the  k-factors  referenced  to  the  dynamic  pressure  remain  relatively 
constant  for  many  restriction  components  over  a  wide  range  of  flow 
conditions.  Therefore,  when  information  does  not  exist  to  permit  the 
extension  of  a  model  to  a  broad  spectrum  of  operating  environments,  the 
application  of  the  incompressible  characteristics  to  high  velocity  flows  still 
may  be  warranted . 
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II.  THEORETICAL  ANALYSIS 


The  analyses  of  the  internal  flow  systems  of  gas  turbine  engines  are  based  on 
traditional  one-dimensional  compressible  formulas  where  the  parameter 
gradients  and  the  velocity  profiles  are  approximated  by  "average"  conditions, 

V  =  m  /  p  A 

Then  the  apparent  loss  of  "average"  total  pressure,  which  results  from  the 
use  of  effective  velocity  to  represent  velocity  profile,  is  absorbed  in  the 
real  system  total  pressure  loss  (1).^  The  internal  flow  system  is  "complex" 
from  the  standpoint  that  most  of  its  duct  geometry  changes  (form  drag)  tend  to 
overwhelm  wall  friction  (skin  drag)  as  the  source  of  total  pressure  losses. 

The  proximity  of  the  restrictions  in  series  is  such  that  fully  developed 
laminar  or  turbulent  flow  is  seldom  achieved.  This  transitional  flow 
environment  contributes  to  the  uncertainty  of  the  one— dimensional  analysis. 
However,  the  basic  simplicity  of  the  formulation  and  the  ability  to  iterate 
the  model  coefficients  from  experience  make  the  approach  viable  and  perhaps 
preferable . 

The  calculation  for  internal  flow  system  performance  generates  a  network  of 
flows,  pressures,  and  temperatures  throughout  the  cooling  and  leakage  paths  in 
the  engine  The  steady-state  solution  has  an  electrical  analog  where  tne 
restriction  k-factors  are  similar  to  the  resistances  and  the  flows  (currents) 
are  found  by  Kirchoff's  law.  The  total  pressures  are  comparable  to  the  node 
voltages.  The  most  important  descriptors  for  a  typical  internal  flow  system 
model  are  the  accurate  values  for  flow  path  cross-sectional  areas  The  flow 
area  is  crucial  to  the  determination  of  flow  and  is  of  first  order  importance 
in  the  estimation  of  the  total  pressure  loss  across  restrictions .  The  basis 
for  defining  the  k-factor  models  and  the  application  of  these  models  to  the 
duct  geometries  which  are  encountered  in  gas  turbine  internal  flow  systems  are 
discussed . 

3  Bracketed  numbers  refer  to  References,  page  69.. 
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Selection  of  the  Generalizing  Parameter  for  k-Factor 


Conventional  usage  employs  the  dynamic  pressure, 
q  =  pV2  /  2  , 

as  the  reference  parameter  for  generalizing  internal  as  well  as  external  drag 
and  pressure  loss  coefficients  for  incompressible  flows.  However,  when 
compressibility  effects  become  important  to  high  velocity  flows  of  gases, 
there  seems  to  be  no  consensus  on  the  reference  parameter  for  drag  and 
pressure  loss  calculations.  External  aerodynamics  has  retained  the  dynamic 
pressure  reference  on  a  uniform  basis, 

Fd  =  Cd  <Vo/2)Aref 

Internal  aerodynamics  vacillates  between  the  traditional  dynamic  pressure,  q 
and  the  impact  pressure,  (P  -  p) , 


P_  =  k  q 
2  max 


P1  -  P2  =  k  <P  -  P>max 


Both  parameters  are  functions  of  Mach  number  and  Y, 


-a-  =  JLn* 

P  2 


(  1  + 


Y  -  1 


Y 

1  -  Y 


Y  -  1  2 

=  1  -  (  1  +  -  MZ  ) 


Y 

1  -  Y 


Therefore,  either  is  capable  of  serving  as  the  generalizing  parameter  for  the 
kinetic  energy  effects  in  compressible  flow.  In  fact,  the  pressure  loss 
coefficient  based  on  reference  area,  A^,  for  any  restriction  can  be 
converted  between  a  dynamic  pressure  basis  and  an  impact  pressure  basis 
without  loss  of  accuracy  or  generality  as 


Y  -  1  2  Y  -  1 

< 1  +  -r  "n  > 


For  incompressible  flow. 


=  k 


since 


P  -  p  =  q 


Some  k-factor  data  are  referred  to  restriction  exit  (downstream:  n  =  2) 
conditions.  This  k-factor  definition  requires  an  additional  iteration  to 
establish  the  total  pressure  loss  across  each  restriction.  For  a  given 
flowrate,  total  temperature,  and  duct  area,  a  downstream  total  pressure  must 
be  assumed. 


When 

P2  ~  P1  k2  q2  =  P2  assumed 


the  solution  has  been  found.  This  complication  is  avoided  by  using  a  direct 
serial  solution  involving  the  maximum  dynamic  pressure  at  the  restriction 
entrance  (upstream:  n  =  1)  as  the  reference  parameter.  The  upstream  and 

downstream  k-factors  for  a  specific  restriction  are  related  implicitly  by 


A  certain  commonality  is  exemplified  by  the  conventional  use  of  dynamic 
pressure  as  the  reference  parameter  for  the  surface  drag  coefficient  for  the 
compressible  flow  over  immersed  bodies  and  for  the  wall  friction  in  conduits 
(FANNO  flow). 


Benedict  and  Carlucci  (2)  have  shown  that  the  application  of  k-factor  values 
based  on  inlet  conditions,  k^,  to  the  equivalent  length  analysis  where 


will  overestimate  the  total  pressure  loss  for  compressible  flow.  1  he  k-factor 
or  4  f  £/  0  in  this  FANNO  type  analysis  is  applied  to  a  continually  increas¬ 
ing  dynamic  pressure  due  to  a  uniformly  distributed  loss  mechanism  through  the 
"constant  area"  restriction.  Therefore,  a  smaller  k-factor  correlates  with 
the  total  pressure  loss  at  a  given  inlet  flow  condition. 

Very  little  reliable  data  are  available  on  the  k-factors  for  restrictions  of 
conventional  geometry  operating  in  the  compressible  flow  regime.  Almost  no 
data  exist  for  the  more  unusual  restriction  configurations  common  to  gas 
turbine  internal  flow  systems.  The  best  data  are  normally  found  for 
incompressible  flow  through  typical  pipe  and  duct  geometries.  Fortunately, 
the  k-factors  based  on  maximum  inlet  dynamic  pressure  are  relatively 
insensitive  to  Mach  number  for  many  restrictions  (3).  The  compressibility 
effects  generally  become  important  above  Mach  0.3  where  the  velocities 
approach  or  exceed  the  critical  Reynolds  number  so  that  the  flow  is  in  the 
fully  turbulent  regime. 

The  selection  of  maximum  inlet  dynamic  pressure  as  the  reference  parameter  for 
total  pressure  loss  coefficients  has  the  merit  of  minimizing  coefficient 
sensitivity  to  compressibility  effects  for  most  loss  mechanisms.  The  dynamic 
pressure  is  analogous  to  the  kinetic  energy  of  the  fluid  stream.  The  impact 
pressure  includes  the  latent  energy  absorbed  by  the  compressibility  of  the 
fluid  in  addition  to  the  kinetic  energy.  The  maximum  inlet  dynamic  pressure 
was  chosen  as  the  reference  parameter  for  generalizing  the  characteri sties  of 
the  total  pressure  loss  coefficients  for  all  of  the  restriction  geometries 
investigated,  with  the  exception  of  the  sudden  expansion.  The  use  of  the 
maximum  inlet  impact  pressure  to  characterize  the  sudden  expansion  results  in 
the  advantage  of  a  unity  loss  coefficient  for  any  jet  discharging  into  a  large 
plenum.  The  selection  of  dynamic  pressure  or  impact  pressure  as  the 
generalizing  parameter,  and  the  choice  of  reference  area  at  the  inlet  or  exit 
of  the  restriction  is  arbitrary .  However,  the  k-factor  must  be  applied  to  the 
value  of  the  generalizing  parameter  for  which  it  was  derived  at  the 
restriction  area  to  which  it  was  referenced  to  produce  correct  predictions  of 
total  pressure  loss. 


Most  pressure  loss  data  available  in  the  public  domain  are  for  circular  or 
rectangular  duct  cross-sections.  Some  data  exist  for  annular  duct  geometries, 
but  their  restriction  configurations  are  limited  primarily  to  constant  area 
and  gradual  expansions.  However,  many  unusual  duct  shapes,  and  particularly 
annular  ducts,  are  encountered  in  the  analysis  of  the  internal  flow  systems  of 
gas  turbine  engines.  Consequently,  when  the  diameter  of  an  equivalent 
circular  cross-section  is  required  for  the  evaluation  of  a  flow  parameter  such 
as  Reynolds  number  or  equivalent  duct  length,  the  hydraulic  diameter  is 
generally  employed. 


HD  = 


4  A 


For  a  circular  duct  HD  =  d 


A  =  -7-HD2 
4 


For  an  elliptical  duct  HD  ~ 


2  a  b 


V2  (  aZ  +  bZ  ) 


so  that 


so  that 


4  2 


For  an  annular  duct  HD  =  D  -  d 


so  that 


(  D  +  d  \ 

A  =  if  \—  j 


HD 

2  /  2 


For  a  rectangular  duct  HD  = 


2  a  b 


C  a  +  b  ) 


so  that 


A  =  (  a  +  b  ) 


HD 


An  analogy  between  annular  and  rectangular  cross-sections  reveals  that  as  a/b  -* 
the  rectangle  becomes  similar  in  geometrical  characteristics  to  an  annulus  of 
small  hydraulic  diameter  where 


/  D  +  d  \ 

\  2  )  a  +  b  ) 


10 


V-'lv!-.' 


* 

I' 

«< 


a/b  -»  0 


or  more  simply 


tr  D  =  b 


and  1 


(  D  -  d  )  =  a 


This  artifice  permits  the  evaluation  of  many  annular  restrictions  from  the 
comparable  restriction  data  for  the  analogous  rectangular  duct. 

The  applications  of  the  theoretical  analyses  selected  for  deriving  pressure 
loss  algorithms  for  turns  and  bends,  combining  and  dividing  branches,  sudden 
expansions  and  contractions  of  flow  area,  and  orifices  will  be  discussed  in 
the  following  sections. 


r/r.H  v 


III.  TOTAL  PRESSURE  LOSS  COEFFICIENTS  FOR  TURNS  AND  BENDS 


Bends  are  among  the  more  difficult  internal  flow  loss  geometries  to  estimate 
accurately  The  duct  geometry  and  condition  of  the  flow  exert  very  strong 
influences  on  the  total  pressure  loss  due  to  the  generation  of  complex 
secondary  flows  and  downstream  recovery  processes.  For  example,  circular-arc 
bends  of  round  or  square  cross-section  develop  twin  counter-rotat 1 ng  helical 
vortices  which  tend  to  stabilize  the  flow,  as  shown  in  Figure  1.  If  the  duct 
cross-section  is  unconventional,  eg.,  triangular,  polygonal,  etc.,  the 
secondary  flow  can  become  complicated  with  more  than  two  vortices. 

Inversely,  turning  in  a  duct  with  a  narrow  annular  cross-section  may  not 
produce  any  secondary  flow.  Combining  the  effects  of  duct  shape  and  wall 
roughness  with  the  rate  and  amount  of  turning  makes  the  loss  analysis  for 
simple,  single  circular-arc  bonds  very  difficult 

The  location  and  flow  environment  of  bends  in  turbine  engine  cooling  and 
leakage  paths  rarely  meet  the  modeling  criteria  for  upstream  and  downstream 
tangent  lengths  or  fully  developed  velocity  profiles.  The  total  pressure  loss 
in  a  bend  is  very  sensitive  to  the  conditions  in  the  entering  flow  as 
established  by  upstream  tangent  length,  wall  roughness,  and  flow 
disturbances  The  length  of  the  downstream  tangent  and  flow  blockage  is  even 
more  important  to  the  pressure  loss  as  the  result  of  the  nature  of  the 
recovery  process  in  the  flow  leaving  the  bend.  Bends  in  gas  turbine  engine 
flow  systems  are  routinely  in  the  region  of  influence  of  upstream  and 
downstream  restrictions ,  which  contributes  to  the  difficulty  of  predicting 
total  pressure  losses  In  addition,  the  flow  area  through  the  bends 
frequently  changes  in  cooling  and  leakage  paths.  The  cons iderat ion  of  these 
application  variables  make  the  total  pressure  loss  prediction  for  internal 
flow  system  bends  uncertain  at  best 


Bends  can  be  classified  according  to  distinct  physical  characteristics  of  the 
flow.  The  rate  of  turning  has  the  greatest  effect  or.  the  flow  through  bends 
and  is  used  in  Table  1  to  distinguish  among  the  primary  types  and  the  physical 
processes  dominating  their-  particular  flow  fields. 

Table  I 

Classification  of  Simple  Circular-Arc  Bends  on  the 
Basis  of  the  Loss  Mechanisms  Dominating  the  Flow  F lelds 


Bend 

Type 

r/h 

Predominant  Loss  Mechanism 

Long 

Bend  s 

>  14 

Wall  friction 

Short 

Bends 

<  14 

Combined  flow  separation 

>  0  .  S 

and  wall  friction 

Mitre 

Bends 

<  O.b 

Flow  separation 

Reference  to  Figure  2  shows  that  turning  flow  in  most  gas  turbine  engine 

restrictions  resides  in  the  short  bend  (and  mitre  bend)  category.  As  the 

result  of  the  generally  elevated  pressures  and  temperatures  and  the  high  flow 

velocities  in  gas  turbine  engines,  the  flow  in  short  bends  will  usually  be 

turbulent,  R  >  20000,  as  predicted  by  Figure  3. 
e 

The  analysis  of  internal  flow  systems  in  gas  turbine  engines  must  attempt  to 
account  for  the  effects  of  the  many  variables  which  influence  the  total 
pressure  loss  in  a  bend.  The  effects  of  the  following  parameters  on  bend 
losses  are  usually  considered  where  the  availability  of  quantitative  data 
permi t : 


Bend  Geometry 

cross-sectional  shape 
turning  rate 
amount  of  turning 
area  change 

Bend  Flow  Conditions 

laminar,  transitional,  or  turbulent 
wall  roughness 
upstream  tangent  length 
downstream  tangent  length 


FIGURE  3.  TRANSITION  REGION  FOR  FLOW  IN  LONG  CIRCULAR-ARC  BINDS 
Reference  (4) 


The  evaluation  of  these  effects  will  be  discussed  as  they  apply  to  the  single 
circular-arc  bend  types  The  modeling  of  compound  bends  and  multiple  bends 
will  not  be  considered.  The  bibliography  contains  several  references  which 
treat  these  subjects  in  varying  degrees  of  analytical  depth  Miller  (b)  and 
ESDU  (6)  are  recommended  sources  of  reliable  performance  data  The  algorithm' 
presented  for  bend  losses  apply  to  turbulent,  incompressible  flow  The 
reliable  data  references  seem  to  agree  upon  the  application  of  the 
incompressible  loss  factors  to  compressible  flows  within  the  present 
state-of-the-art .  The  experimental  and  analytical  data  are  not  yet 
sufficiently  reliable  to  warrant  a  distinction  at  this  time 

When  a  more  precise  solution  is  required,  and  the  bend  installation  and  flow 
quality  justify  the  analytical  complexity,  reference  (4)  can  be  used  to  model 
bend  K-factors 

Basic  Circular-Arc  Bends 

The  empirical  model  for  K-factors  proposed  by  Ito  (7)  is  recommended  for 
turbulent  flow  through  bends  of  circular  cross-section  in  references  (1),  (4), 

and  (8),  among  others  The  total  pressure  loss  predicted  by  the  Ito  model  is 
slightly  greater  than  that  predicted  in  reference  (9).  However,  the  wall 
friction  loss  is  included  in  the  Ito  formulas  while  the  plots  in  reference  (9) 
represent  turning  loss  alone, 

Kb(9)  "  Kb(7)  "  Af  6  (~HD  ^ 

The  bend  model  by  Ito  is  limited  to  hydraulically  smooth  walls  so  that  a 

correction  for  rough  walls  is  required.  Although  the  algorithm  has  been 

4  b 

validated  by  test  data  for  a  Reynolds  number  range  of  2  (10  )  to  4  (10  ), 

the  formulas  can  be  extrapolated  to  a  Reynolds  number  of  1  (10^)  with 

acceptable  accuracy  Bend  loss  does  not  change  significantly  with  Reynolds 

6 

number  greater  than  1  (10  ) 


Long  Bonds  Kg  ^  ^  <  3b4 

K.  -  0  01746  a  f  0  (-F-  ) 
b  c  \  h  / 

Since  this  equation  applies  to  a  minimum  r/h  of  7.4,  its  use  in  gas  turbine 
internal  flow  analysis  arises  infrequently. 


The  secondary  flows  present  in  bends  of  circular  cross-section  generate  addi¬ 
tional  losses  due  to  wall  friction.  These  smooth— wall  friction  losses  are 

correlated  for  curved  turbulent  flow  at  Re  CO  <  1200  by  H  Ito  (10)  as 


f 

c 


0.0205  0. 304 


/' 


The  implied  region  of  validity  extends  to  short  bends  with  r/h  as  small  as  41 
in  addition  to  the  entire  long  bend  envelope. 


Short  Bends  (not  including  mitre  bends) 


-o.  i; 


K  =  0.00431  a  e  Re 


(v) 


0.84 


The  turning  losses  for  long  and  short  bends  are  correlated  by  the  a  term  as 
determined  from  Table  II. 


A  linear  interpolation  between  the  defined  bend  angles  produces  consistently 
smooth  k-factor  charac ter l st i c s ,  as  illustrated  in  Figure  4.  The  bend  loss 
charac ter i s t l c s  generated  for  90  deg  short  bends  are  plotted  in  Figure  5. 


Idel'chik  contends  in  reference  (13)  that  all  turns  and  bends  are  essentially 
independent  of  the  relative  roughness,  e/  D,  of  the  wall  at  Reynolds  numbers 


Table  II. 

Turning  Loss  Factors  for  the  Bend  Loss 
Model  by  I to  ( 7 ) , 
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FIGURE  5  SINGLE  SHORT  C 


less  than  4  (10  ).  The  formulation  for  the  k-factor  of  long  bends  suggests 
a  direct  ratio  of  friction  factors  for  a  wall  roughness  correction  as 


_  c  rough 
f  ~  f 

c  smooth 

Unfortunately,  reliable  data  to  evaluate  or  substantiate  this  hypothesis  were 
not  found  for  friction  factors  of  curved  flows.  Miller  (5)  suggests  such  a 
wall  roughness  correction  based  on  straight  pipe  friction  factors  applied  to 
all  circular-arc  bends.  Idel'chik  (13)  restricts  this  correction  factor  to 
circular-arc  bends  with  r/ HD  <  1.5.  Both  Idel'chik  (13)  and  Henry  (14) 

5 

propose  a  stronger  influence  of  wall  roughness  at  Reynolds  numbers  above  2  (10  ) 
and  for  short  bends  with  rt HO  >  1.5.  However,  a  maximum  effect  of  C^.  -  2.0 
is  proposed  by  Idel'chik  for  any  combination  of  wall  roughness  and  Reynolds 
number.  An  approximate  model  combining  this  general  concensus  was  synthesized 
for  short  circular-arc  bends  as  shown  in  Table  III. 


Table  III. 

Effect  of  klall  Roughness  on  Short  Circular-Arc  Bends  . 


R  <  4  ( 10) A 
e 


Cf  =  1.0 


4  (10A)  <  Rp  <  2  ( 105) 
r/HD  <  1.5 

r/HO  >1.5 
Re  >  2  (105) 


C,  =2.0 
f  max 


C„  = 


rough 


f  f 


Mitre  Bends 


The  circular-arc  bend  degenerates  into  a  special  case  where  the  concentric 
inner  bend  radius  goes  to  zero  at  r/h  =  0.5.  Geometrical  interfaces,  size 
limitations,  or  ease  of  fabrication  produce  many  bend  restrictions  with  corner 
points  at  the  inside  wall  and  outside  wall,  r/h  =  0.  These  bends  and  certain 
variations  of  similar  type  are  categorized  as  mitre  bends.  The  unifying 
characteristic  of  the  flow  through  mitre  bends  is  the  high  rate  of  turning 
The  separation  and  turbulent  mixing  flow  processes  dominate  the  total  pressure 
losses  in  mitre  bends  so  that  Reynolds  number  effects  are  small  to  quite  low 
values.  Some  dispersion  is  noted  from  source  to  source,  but  an  average  of 
data  from  references  (1),  (4),  (9),  and  (11),  plotted  in  Figure  6,  is  a  good 
representat ion  of  the  group.  For  most  internal  flow  systems  found  in  gas 
turbine  engines  the  curve  of  Figure  6  can  be  adequately  reproduced  by  the 
equation  proposed  by  Hager  (12)  for  bend  angles  greater  than  25  degrees, 

\  =  2  (1  -  COS 

As  a  rule  of  thumb,  between  bend  angles  of  5  degrees  and  25  degrees 


K 


b 


K 

b(  12) 


+  0.02 


can  be  used . 


Figure  7  shows  the  evolution  of  the  mitre  bend  and  some  of  the  variations 
which  are  encountered  in  practice,  Experience  with  these  modified  mitre  bends 
has  shown  that  the  radius  on  the  inside  wall  is  the  most  influential  geometry 
factor  for  reducing  the  K-factor  (total  pressure  loss). 

Corrections  applicable  to  Figure  6  for  these  geometric  variations  are  provided 
in  F igure  8 . 
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PRESSURE  LOSS  FACTOR 


SHORT  CIRCULAR 
r/h  0.5 


DEGENERATE 

CIRCULAR-ARC 

r/h  =  0.5 


I 


h 


MITRE 


RADIUS,  R;  INNER  RALE  RADIUS 
R0  OUTER  WALL  RADIUS 

FIGURE  8  TOTAl  PRESSURE  I  OSS  CIIARAf  TfRfST fCS  FOR  nopiFIED  MIJRF 
OF  CONS  I  AN  I  AREA  Willi  CUNSIANT  HEIGHT,  ,  AT  INLET  AND 


The  data  of  references  (A)  and  (13)  indicate  that  very  rough  walls  car, 
increase  the  k-factor  for  mitre  bends  as  much  as  SOX  at  Reynolds  numbers  above 

4(104), 

cf  =  1  ♦  5  “°3>(hs) 

where  C  ~  1 . 5 

f  max 

4 

Below  a  Reynolds  number  of  4  (10  )  the  effect  of  roughness  or  Reynolds 
number  on  mitre  bend  k-factor  is  negligible. 

Effects  of  Other  Geometrical  Parameters 

The  basic  bend  model  was  derived  for  ducts  of  circular  cross-section,  but 
testing  has  shown  that  the  algorithm  represents  flow  through  square  ducts 
almost  as  well  From  this  point  the  bend  model  can  be  extended  to  include 
ducts  of  elliptical  and  rectangular  cross-section.  The  bend  model  was 
formulated  using  the  geometrical  parameters  0  and  h  to  accommodate  this 
extended  scope. 

The  pressure  drop  factors  for  the  flow  through  bends  of  square  and  circular 

cross-section  at  the  same  values  of  r/h  (including  mitre  bends),  0,  and  R^ 

are  assumed  to  be  negligibly  different.  Experimental  correlations  by  Ward-Smith  (IS) 

S  S 

and  Hiller  (16)  demonstrated  this  premise  for  1  (10  )  <  Rp  <  13  (10  ). 

Figure  9  is  a  correlation  of  rectangular  duct  performance  relative  to  circular 
ducts  due  to  Ower  and  Pankhurst  (17).  The  correction  for  rectangular  bonds 
can  be  represented  by 
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Rcuunce  of  beods  in  pipes  of  rcosnguUr  cross-section. 


k  -  FACTOR  FOR  TOTAL 
PRESSURE  LOSS  IN  A 
BEND  OF  SQUARE  CROSS- 
SECTION,  b  =  h. 

k  =  FACTOR  FOR  TOTAL 
PRESSURE  LOSS  IN  A 
BEND  OF  RECTANGULAR 
CORSS-SECTION  WITH 
ASOECT  RATIO  b/h. 

FIGURE  9  RESISTANCE  COEFFICIENT  FOR  BENDS  OF 

RECTANGULAR  CROSS-SECTION.  REFERENCE  (17) 


2b 


where 


b/h  <1  >1 

m  0.2551  0.1386 


The  performance  of  elliptical  ducts  can  be  estimated  from  this  correlation  as 
a  qualitative  approximation.  If  more  accurate  analysis  is  required,  reference 
(9)  or  (14)  can  be  consulted  for  extensive  data  on  bend  performance  of 
elliptical  and  rectangular  ducts. 


The  basic  bend  model  was  derived  from  tost  data  for  fully  developed  turbulent 

flow  at  the  inlet  to  the  bend  and  for  a  least  fifty  diameters  of  downstream 

tangent  length.  The  downstream  tangent  generally  contributes  mixing  losses 

for  lengths  greater  than  two  diameters.  However,  as  the  outlet  tangent  length 

diminishes  toward  zero,  the  k-factor  increases  as  the  initial  pressure 

recovery  process  in  the  first  two  to  four  diameters  of  downstream  tangent  is 

lost.  Miller  (5)  provides  a  convenient  correction  for  downstream  tangent 

length,  shown  in  Figure  10.  A  k*  =  1.2  curve  based  on  data  in  reference 

b 

(4)  has  been  added  ostensibly  for  corrections  to  mitre  bend  k-factors .  For 
ducts  of  particular  rectangular  cross-section  Miller  recommends  the  following 
modifications  to  of  Figure  10: 

if  b/h  <  0.7  and  fi.  /HD  >  1, 
d 


1  +  C 


~fi.r 


If  the  bend  or  downstream  tangent  discharges  into  a  larger  duct  or  plenum,  a 


sudden  expansion  loss  must  be  added  to  C 


Note  that  neither  the  bend 


k-factor  nor  include  the  wall  friction  loss  (4f  £./  D)^  associated 


with  the  downstream  tangent  length. 


The  flow  in  gas  turbine  internal  flow  systems  is  often  turned  in  annular 

ducts.  These  annular  flow  paths  and  conventional  duct  paths  incorporate 

turning  with  area  change  in  many  restrictions  The  effect  of  turning  with 

area  change  in  circular  arc  bends  was  correlated  by  Henry  (14),  Figure  11 

Data  for  k-factors  of  mitre  bends  with  area  change  were  compiled  by  Idel'chik 

(13).  A  similar  presentation  of  the  data  is  provided  in  Figure  12  The 

severe  effect  of  the  sharp  corner  on  the  inside  wall  is  evident  from  the 

inversion  experienced  by  the  loss  coefficient,  C  .  A  sudden  expansion 

ft 

alleviates  some  of  the  restriction  at  the  mitre  corne1'.  The  flow  turns  more 
gradually  with  less  loss  A  contraction  at  the  separation  point  of  the  mitre 
corner  restricts  the  downstream  area  available  for  turning.  The  flow 
contraction  is  amplified  and  the  k-factor  increases. 


WIDTH  OF  ENTRANCE/WIDTH  OF  EXIT,  yh2 


FIGURE  11  TOTAL  PRESSURE  LOSS  COEFFICIENT  FOR  SO  DEC 
CIRCULAR-ARC  BENDS  OF  CHANGING  AREA, 
REFERENCE  («) 


TOTAL  PRESSURE  LOSS  COEFFICIENT 


IV.  TOTAL  PRESSURE  LOSS  COEFFICIENTS  FOR  BRANCHES 


The  steady  flows  through  junctions  and  branches  differ  from  flows  through 
other  restrictions  discussed  in  that  mass  is  increased  or  decreased  within  the 
component.  Most  restrictions  constitute  series  flow  losses  where  the  massflow 
leaning  is  the  same  as  the  massflow  entering  However,  internal  flow  systems 
contain  many  flow  intersections  of  parallel  restrictions  where  the  local 
dynamics  influence  the  total  pressure  loss  The  flow  model  for  the 
intersection  should  include  k-factors  which  account  for  the  effects  of  mixing 
in  combining  flows  and  diffusion  turbulence  in  dividing  flows  These  flow 
processes  are  generally  a  function  of  the  flow  split,  which  requires  an 
iterative  solution  for  the  correct  k-factor.  The  flow  processes  in  junctions 
and  branches  behave  like  those  in  bends  in  many  ways,  but  the  additional 
effects  of  mixing  flows  from  different  sources  or  delivering  flows  to 
different  sinks  complicate  their  physical  models  Multitudes  of  junction  and 
branch  geometries  exist  in  engineering  practice.  Two  of  the  most  common 
geometries,  symmetrical  and  unsymmetrical ,  are  shown  in  Figure  13 

fit  a  junction  or  a  branch,  continuity  must  be  satisfied 


m3  =  m2  +  mj 


In  addition  momentum  and  energy  must  be  conserved 


P3  A3  4  P3  A3  V3  ~  ^P1  fll  4  P1  P1  'V  C°5  6  4  ^P2  A2  4  P2  °2  V2^  C°S  P 


where 


cos  6=0  for  the  unsymmetrical  case 


me  T=mc  T+mc  T 
3  P3  3  1  Pj  1  2  P2  2 


(ad labat ic ) 


34 


• .  •  .  ■  ,  >  ,  •  .  v  v  >  ■/  v-v;  •  /.  r.  , 


. 


,  . 


Sft 


> 


FIGURE  13  COMMON  GEOMETRIES  FOR  JUNCTIONS  AMD  BRANCHES 
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Vazsonyi  (18)  attempted  the  analytical  prediction  of  junction  and  branch 
performance  on  the  basis  of  this  physical  model  utilizing  bend  flow  process 
analogies.  The  adaptation  of  this  work  for  generalized  junction  and  branch 
k-factors  by  the  SAE  (9)  was  unsuccessful  Although  a  profusion  of  pressure 
loss  data  exist  for  a  multitude  of  geometries,  the  generalization  and 
reduction  to  usable  k-factor  parametrics  is  limited  essentially  to  the  basic 
symmetrical  and  unsymmetrical  varities.  Williamson  and  Rhone  (19),  however, 
do  present  a  survey  of  special  geometries. 


The  k-factors,  k  ,  have  been  defined  using  the  dynamic  pressure  in  the 
n  3 

limb  with  the  combined  flow  (leg  3).  The  k-factor  is  positive  for  a  total 
pressure  loss  or  is  negative  for  a  total  pressure  gain.  The  k-factors 
presented  for  junctions  and  branches  do  not  include  the  total  pressure  lost  to 
wall  friction.  Consequently,  the  loss  model  for  combining  or  dividing  flows 
with  upstream  tangent  and  downstream  tangent  at  least  three  hydraulic 
diameters  long  should  be  as  follows: 


n°3  P  =  Vl)+  (  2  A)*  4fn(ii)  (  2  Pn  “n  ) 


Occasionally  it  is  desirable  to  reference  the  k-factor  to  another  leg  (n), 

r  T2 


As  for  bends,  experimental  results  for  circular  ducts  and  square  ducts  show 
negligible  difference  (1).  Little  influence  from  Reynolds  number  is  evident 
in  turbulent  flow.  When  the  flow  is  not  turbulent,  the  energy  contribution  is 
generally  minimal  due  to  small  dynamic  pressure. 


Although  the  following  junction  and  branch  models  are  based  on  experiments 
with  fluids  at  constant  densities,  the  results  can  be  applied  to  compressible 
flows  with  reasonable  accuracy. 


If  more  analytical  precision  is  required  and  the  restriction  geometry  and  flow 
environment  warrant,  reference  (20)  for  combining  flow  junctions  and  reference 
(21)  for  dividing  flow  branches  can  be  used  for  restriction  modeling 


Symmetrical  Junctions  and  Branches 


Symmetrical  junctions  and  branches  are  often  referred  to  as  wyes  because  of 
their  geometrical  configuration.  The  data  of  Miller  (16)  correlate  well  with 
that  of  other  investigators  and  are  slightly  pessimistic.  It  is  ordinarily 
good  design  practice  to  err  on  the  side  of  high  total  pressure  loss,  so  the 
Miller  (5)  performance  maps  in  Figure  14  were  selected  for  the  total  pressure 
loss  associated  with  combining  flows  in  wyes. 


lk3 


P1  - 


2  P3 


P 


V 


2 

2 

3 


Simi larly 
the  total 


,  the  performance  maps  by  Mi  lie 

pressure  loss  model  of  dividing 

P  -  P 
3  1 


2  P3  V3 


(5)  in  Figure  15  were 
flows  in  wyes, 


pref erred 


for 


Another  common  class  of  symmetrical  junction  is  the  90°  three-way  dividing 
branch  or  four-way  cross.  Miller  (5)  provides  k-factor  maps,  shown  in  Figure 
16,  for  the  perpend l cu lar  off-take  leg,  jkg.  anc*  for  the  straight-through 
leg,  2^3'  The  performance  is  shown  for  a  dividing  junction  with  all  legs 
of  equal  area  and  with  sharp  edges  at  the  intersections. 


Unsymmetrical  Junctions  and  Branches 

Junctions  and  branches  having  two  of  the  limbs  colinear  are  frequently 
encountered  in  gas  turbine  internal  flow  systems.  Restrictions  interfacing  at 
90  deg  tees  are  common.  The  modeling  of  manifolds  is  one  of  the  most 
important  applications  for  such  k-factor  data.  Fortunately,  Gardel  (22)  has 
done  a  comprehensive  experimental  program  to  determine  the  effects  of 
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FIGURE  15  DIVIDING  FLOW- SYMMETRICAL  T  JUNCTION 
REFERENCE  (5) 


Ion  coefficient 


changing  duct  sizes,  lateral  angles,  and  introducing  fillets  and  radii  at  the 
lateral  limb  Gardel  derived  the  following  empirical  equations  to  represent 
the  performance  of  unsymmetrical  junctions  (combining  flow): 


.2  _  2 


lk3  -  0.92  (i  -jr  [(i.2  -r  X 


1/2  w  cos  0 


+  0  8  ( 1 


-  jg)  -  (^-  1)  cos  b]+  (2  -a.)(  1  -  J)  <j£ 


2k3  =  0.03  (1 


.  .  -  038  (1 


-ef)2  -<j^  [i  +  (i.62  -r1/2)( 

a.)]  *  ( 2  -  du  a  ~f) 


1/2 . .cos  6 


where 


P  -  P  P  -  P 

13  2  3 

k  =  -1 -  and  k  =  - 2 

1/2  P3  V3  1/2  p3  V3 


for  unsymmetrical  branches  (dividing  flow): 


k3  =  0.95  (1  [(1. 


2  +  [(1.3  tan  |  -  0 . 3  4 


O'0-* ‘£)/2)J,0-4t<I-f 


2k3  =  0.03  (1  -^J)2  4  0  35^T  -  02^(1  -^) 


where 


P3  -  P1  P3  P2 

k  =  — -  -  and  k  =  — - - — 

1  3  2  2  3  2 

1/2  P3  V3  1/2  p3  U3 


PVi,Vi 


■■  ‘  .•/  ^V.'  v’v’ 

..  •  m  “>  a  ,*  a  a"*  AjiN  hf  m  k *  «  *  a  “  a  *  •  *  k  « 

A  VL- .M  w  iLlk.  JL-.  Km  mV.  A  -  A*  *  a  _  a*_ 


In  these  equations 


for  a  range  of 

lateral  angle 

A 

CL  r 

ib°  <  e 

Also 

(0  62b 

and 

(0  <  | p 

<  165° 

A»<  1) 

0.12) 


For  compressible  flows  Qf= 


is  preferred 


V.  TOTAL  PRESSURE  LOSS  COEFFICIENTS  FOR  SUDDEN  AREA  CHANGES 


The  most  common  restrictions  encountered  in  modeling  internal  flow  systems  for 
gas  turbine  engines  are  sudden  expansions  and  sudden  contractions  as 
illustrated  in  Figure  17. 

The  sudden  expansion  loss  is  well  represented  by  a  one-dimensional  analysis. 
Although  the  sudden  contraction  appears  to  be  the  geometrical  reverse  of  the 
sudden  expansion,  it  is  not  possible  to  obtain  a  comparable  explicit  solution 
for  total  pressure  loss  from  a  one-dimensional  flow  model. 

Sudden  Expansion 

Flow  from  a  duct  into  a  sudden  enlargement  can  be  analyzed  by  conserving  mass, 
momentum,  and  energy  between  the  discharge  plane,  A^  ,  and  the  reattachment 
plane,  A^ .  Incompressible  turbulent  flow  is  well  represented  by  the 
Borda-Carnot  relation 


where 

1°2  p  =  kse  Pi 

Laminar  flow  correlates  poorly  with  the  one-dimensional  equations  due  to. 

o  the  large  velocity  gradients  in  the  profile  of  the  efflux, 

o  the  important  shear  stress  contribution  to  reattachment  with  the  possible 
unsymmetrical  flow  fields  in  two-dimensional  duct  geometries. 
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SUDDEN  EXPANSION 


SUDDEN  CONTRACTION 


FIGURE  17  GENERAL  CONFIGURATION  OF  SUDDEN  AREA  CHANGES 


Adiabatic  compressible  flow  can  be  modeled  similarly  as  described  by  Benedict 
et .  al . ,  (23) 


Y_i  2  1/2  Y-i  2  1/2 

M  (1  +  —  VC)  M  (1  +  M  ) 

2  2  2  1  _ 2  1 


1  +  Y 


1  ♦ v  "i  +(r)(LtA) 


where 


p  -  static  pressure  at  the  face  of  the  step, 
e 


When  the  efflux  from  duct  A  is  subsonic  (0  <  M  <  1),  the  flow  field  in 

1  ~  1  - 

the  enlargement  is  also  subsonic  (0  <  M  <  1)  and  p  /p  -  1.0. 

~  2  ~  e  l 


Ward-Smith  (1)  suggests  using  a  parameter 


Y-l  2  ' 
(1  +  ~  M  ) 


which  simplifies  the  subsonic  equation  to  the  recognizable  quadratic  form 
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The  total  pressure  ratio  is  determined  by  the  expansion  geometry  and  Mach 
numbers  as 


*  ■  tl) 


Y-l  2 

1  +  T  M2 


Y-l  2 
1  +  —  M1 


y+  i 


2  (Y-l) 


Then  either  the  k-factor  based  on  dynamic  pressure 
P. 


1  "  P1 


Y  / ^l\ 

2  1  P  I 


or  the  k-factor  based  on  impact  pressure 
P 


se 


1 


1  -  r 

i 


can  be  calculated. 


Singularities  are  encountered  where  becomes  very  large  (A.  -+  0)  and 


whore  M  becomes  very  small  (M  -»  0) .  Noting  that  N  and  M  become 
1  1  a  2  2 


zero  when  \  -  0  for  any  value  of  M  the  entire  energy  in  the  jet,  (P^ 


is  dissipated  in  the  expansion  to  p£  =  pj .  Since  k$=  1.0  at  A  =  0, 


1  -  (Pj/Pj) 


SP 


5  (■)  • 


at  A,  =  0. 


Numerical  realities  make  the  use  of  k  or  k  =  1.0  at  A 

se  se 


0  recommended 


practice  for  all  A  <  0.0001 


Any  compressible  fluid  flowing  adiabatical ly  assumes  a  constant  density  charac¬ 
ter  at  Mach  numbers  below  0.1.  Consequently,  the  incompressible  Carnot-Borda 
equation  for  a  sudden  expansion  loss  can  be  used  as  the  asymptotic  value  for 


k  se  and  k  se  at  Mi  <  0.1. 
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The  relationship  between  the  k-factors  based  on  dynamic  pressure  or  impact 
pressure  are  shown  in  Figure  18  for  all  subsonic  sudden  expansions  of  a 
perfect  gas  with  Y  =  1.4.  Although  either  k-factor  definition  is  equally 
accurate,  physical  conceptual ization  of  the  loss  as  an  extension  of  the 

incompressible  case  seems  easier  based  on  impact  pressure,  k^  . 

When  the  efflux  from  duct  A^  is  supercrit ical  (M  >  1.0)  for  the  choked 
condition,  PI  =  1.0,  the  sudden  enlargement  equation  can  be  solv/ed  for  the 
effective  step  face-to-jet  static  pressure  ratio  as  a  function  of  PI 


\  yY+lN1 
Pj  ~  \1-^ A  2  ) 


1  +  Y 


2  ,  /  Y+l  \ 

/  Izl  2\1/2  "  \  2  / 
M2  [l  +  2  M2) 


The  total  pressure  loss  calculation  in  this  underexpanded  flow  regime 
simplifies  to 


Y-l  2~1  — 

1  +  —  M2  2  (Y— 1 ) 


n±n 

2 


An  iterative  solution  can  then  be  performed  to  determine  the  PL ,  p  /p 

Z  el 

pair  compatible  with  the  given  duct  discharge  conditions  P  ,  p  .  Assuming 

1  1 

a  final  subsonic  Plach  number,  P I  <  1,0,  the  terminal  conditions  must  conform 


T  \fmr.  h  y+l 
*2  v  R  (,  ♦ 


The  k-factors  for  the  supercritical  sudden  expansion  can  be  easily  determined 
from  the  constant  pressure  ratio  of  the  jet  as 


_ h. 

0  36980 


Since  the  Borda-Carnot  derivation  is  unconcerned  with  the  duct  shape,  the 
quality  of  fit  for  the  sudden  expansion  model  indicates  that  duct  shape  is  a 
secondary  effect,  as  shown  in  Table  IV.  Variations  in  the  discharge  port 
geometry  have  a  negligible  effect  on  the  sudden  expansion  loss.  Even 
gradually  expanding  ducts  with  divergence  angles  exceeding  about  45° 
experience  flow  separation  which  behaves  essentially  like  an  abrupt 
enlargement  The  concept  can  be  extended  to  include  the  separated  flow  (vena 
contracta)  at  sharp-edged  entrances  or  forward  facing  steps.  The  physics  of 
the  sudden  expansion  model  are  helpful  in  understanding  the  modeling  of  the 


k-factors  for  sudden  contractions  and  sharp-edged  orifices. 

Table  IV. 

Less  k-factor  for  pipe  exits.  Reference  (24) 

F  1 1 1 1 nq 

Description 

k-factor 

— - . 

j 

Projecting 

1  .0 

I 

Sharp  edged 

1.0 

Rounded 

1.0 

j 

Sudden  Contraction 

Many  investigators  have  studied  the  flow  modeling  of  a  duct  entrance  from  a 
fluid  reservoir  or  an  abrupt  area  reduction  within  a  duct.  Adiabatic  flow 
from  an  infinite  reservoir  into  a  re-entrant  duct  (Borda  mouthpiece),  Figure 
19,  can  be  accurately  modeled  one— d imens ional ly  by  conserving  mass,  momentum, 
and  energy  between  the  fluid  reservoir  and  the  vena  contracta. 


y 


c 


where 


(*  * 


V-l 


se 


1  - 


P 

_2 

P 

1 


0.47172 


lim  C 


1 


or 


The  total  pressure  loss  in  this  type  of  entering  (acce lerat ing )  flow  is 
relatively  small  because  the  conversion  of  pressure  to  velocity  (as  in  a 
nozzle)  is  a  stable  process.  Consequently,  the  assumption  of  inviscid  flow 
between  the  reservoir  (station  1)  and  the  vena  contracta  (station  c)  is 
appropriate ,  Then  the  sudden  expansion  model  for  compressible  adiabatic  flow 
can  be  employed  to  evaluate  the  total  pressure  loss  for  flow  reattachment 
This  modified  Hughes  and  Safford  analysis  predicts  a  1.0  <  k  <  1.09b  ranqe 

SC 

for  compressible  flow  into  a  Borda  mouthpiece.  Dodge  (24)  reports  a  range  of 
k-factors  from  0.68  to  2.5  for  incompressible  flow  depending  upon  inlet  edge 
conditions  (modified  corners  to  sharp).  This  large  discrepancy  with  empirical 
k-factor  values  suggests  a  shortcoming  of  the  one-dimensional  analysis  as  it 
applies  to  sudden  contraction  losses.  The  marked  curvature  of  the  vena 
contracta  flow  compromises  the  accuracy  of  the  one-dimensional  assumption 


As  the  re-entrant  length  of  the  tube  decreases  to  zero,  the  Borda  mouthpiece 
is  transformed  into  a  sharp— edged  inlet.  The  entrance  flow  no  longer  develops 
in  isolation  from  the  reservoir  wall,  so  the  momentum  analysis  must  be 
modified.  Reaction  of  the  flow  with  the  wall  assists  turning,  and  the 
incompressible  contraction  coefficient  increases  to  about  0.6.  Miller  (5) 
provides  dramatic  data  for  this  effect  as  tube  wall  thickness  increases  for  a 
re-entrant  inlet.  The  momentum  equation  no  longer  defines  the  total  pressure 
loss  explicitly.  The  Hughes  and  Safford  equation 


sc 


2  2 
C  C 
v  c 


C 

c 


predicts  k  =  1.00  for  a  Borda  mouthpiece  and  k 

sc  sc 

edged  inlet  from  the  flow  characteristics  in  Table 


=0.56  for  a  sharp- 
V. 


Table  V. 

Characteristics  for  incompressible  flow  in  duct 
entrances  and  exits,  Reference  ( 25 ) 


Restriction 


C 

v 


C 

c 


C 

D 


re-entrant  inlet  0.98 
sharp-edged  entrance  0.80 
duct  discharge  1.00 


0.52  0.51 
1.00  0.80 
1.00  1.00 
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Benedict,  et .  al.,  (26)  propose  a  generalized  equation  based  on  discharge 
coefficient  to  more  accurately  represent  sudden  contraction  losses  in  constant 
density  flows  with  an  approach  velocity 


-1  1  - 


If  A  /A  is  taken  to  represent  C  for  a  re-entrant  inlet,  k  =  2  08 
2  1  c  sc 

is  predicted,  which  better  represents  the  empirical  value.  A  k  =  0.56  is 

sc 

still  predicted  for  a  sharp-edged  entrance  where  A^  is  very  large  with 

respect  to  the  duct  area.  If  a  baseline  total  pressure  loss  coefficient  is 

defined  for  a  sudden  contraction  as 

*  1 

k  =  —  -  1 

“  i 

then  the  influence  of  the  contraction  ratio  can  be  formulated  as 


k  =  k 
sc  sc 


-fry 


In  many  composite  restrictions  the  flow  at  the  entrance  or  through  a  sudden 
contraction  occurs  at  nearly  constant  density.  Then,  for  most  practical 
applications  to  gas  turbine  internal  flow  systems,  an  incompressible  equation 
is  sufficiently  accurate.  For  A  <  0.3,  which  includes  most  entrances  and 
many  sudden  contractions  within  restr i ct ions ,  the  recommended  equation  is 


k  =  0.5781  (1  -  A  ) 
sc  v  ' 


If  more  accuracy  is  required  or  when  A  >  0  3,  the  least  squares  curve  fit 
from  the  test  data  of  Benedict,  et.  al,  (26)  should  be  used 


k  =  0.57806  +  0  39543  A  -  4.53854  A  ... 

SC  3/2  2  5/2 

...  +  14.24265  A  -  19,22214  A  +  8.54038  A 


These  data  represent  sudden  contraction  characteristics  for  constant  density 

flow  into  long  ducts  where  complete  reattachment  is  assured  If  the  duct 

contraction  length  is  short,  (H/HD)^  <  3,  the  "long  hole"  correction 

presented  under  orifice  restrictions  should  be  applied  to  k  The  effect 

sc 

of  compressibility  increased  the  experimental  value  for  k  as  much  as  12X 
for  subsonic  flow  (26)  Considering  the  uncertainties  associated  with  data, 
installation,  and  environment  the  constant  density  model  for  sudden 
contractions  is  justified 

The  sudden  contraction  model  discussed  so  far  applies  only  to  tubes  or 
entrances  which  are  aligned  with  the  approaching  flow  and  have  sharp  inlet 
edges  figure  20  provides  corrections  obtained  for  entrances  oblique  to  the 
approaching  flow  One  curve  applies  to  a  sudden  contraction  with  the 
downstream  duct  normal  to  the  step  or  wall  which  is  at  an  angle  to  the 
approaching  flow  The  other  curve  applies  to  a  sudden  contraction  with  the 
downstream  duct  at  an  angle  to  the  step  or  wall  which  is  perpend i cu lar  to  the 
approaching  flow.  Figure  21  provides  a  correction  factor  which  accounts  for 
rounding  or  edge  break  effects  on  sudden  contraction  character i s t i c s .  Then 
the  general  contraction  coefficient  can  be  found  as: 

k  =  C  C  k*  (1  -  A2) . 
sc  a  r  sc 

If  better  analytical  precision  is  required,  reference  (27)  for  sudden  duct 
enlargements  and  reference  (28)  for  sudden  duct  contractions  can  be  used  fur 
restriction  modeling 


consist  of  a  thin  end  wa 
between  considerably  lari 


The  flow  charac t er l s t l c  that  distinguishes  an  orifice  from  a  Icing  hole  or  r, 
nozzle  is  the  inability  of  the  vena  contracta,  which  is  induced  by  the 
sharpness  of  the  sudden  contraction,  to  isolate  itself  by  flow  reattachment  to 
the  wall  before  the  sudden  expansion.  Consequently,  the  orifice  flow  model 
must  combine  a  sudden  contraction  and  a  sudden  expansion  which  are  modified  to 
account  for  interactions  with  the  special  flow  processes  that  occur  at  the 
entrance  and  the  exit.  The  conventional  procedure  has  been  to  add  a  modifying 
k-factor,  kj,,  which  incorporates  the  interactive  effects  of  upstream  and 
downstream  duct  geometry  with  the  correction  required  to  account  for  the 
sustained  separation  of  the  flow  through  the  hole.  Consideration  of  the  small 
wall  friction  loss  can  be  included.  Since  the  complex  flow  processes  in  the 
formation  and  disolution  of  the  vena  contracta  are  only  qualitatively 
understood,  the  orifice  model  is  based  upon  empirical  correlations  for 


kn  =  k  +  kn  +  4ffi./H0  +  k 

0  s  C  a.  5  G 

When  the  length  of  the  small  hole  connecting  a  sudden  enlargement  is  less  than 
three  hydraulic  diameters,  the  vena  contracta  formed  at  the  entrance  to  the 
orifice  may  not  reattach  within  the  short  length  Without  reattachment  the 
sudden  contraction  becomes  sensitive  to  flow  conditions  in  the  downstream 
enlargement.  Separated  flow  at  the  orifice  exit  does  not  conform  to  the  model 
established  for  a  sudden  expansion.  An  orifice  model  for  incompressible  flow 
to  account  for  the  process  interactions  caused  by  separtion  was  derived  by 
Dodge  (2A) ; 

$ 

where  k  =0.5 
sc 


* 

where  k^  is  defined  by  Figure  23 


# 

where  k  =1.0 
se 


k  =  k  [ 1— A] 
sc  sc  J 


1/2 


k  =  k  [k  k  ] 
x.  x  sc  se 


k  =  k  [1-71] 
se  se  J 


£/HD 

FIGURE  23  LONG  HOLE  EFFECT.  REFERENCE  (24) 


The  vens  contracts  loss  for  incompressible  flow  through  thick  static  orifices 


is  evaluated  by  the  empirical  k^  of  Figure  23.  7  fie  influence  of  approach 

velocity  and  limited  downstream  expansion  on  the  severity  and  extent  of 
separation  is  corrected  by 

«  1/2 

kJi  =  [1_AJ  [1_^J 

The  performance  of  orifices  is  buffered  against  effects  from  adjacent 

restrictions  or  tangent  geometry  by  1.5  hydraulic  diameters  or  more  of 

straight  duct.  The  orifice  characteristi cs  are  unaffected  by  Reynolds  numbers 

5 

in  the  hole  that  are  greater  than  1  (10  ). 

Fluid  compres s l bi 1 l ty  exerts  the  strongest  influence  on  high  velocity  flow 
through  orifices.  Contrary  to  the  invariant  behavior  of  the  vena  contracta  of 
l ncompre s s i b 1 e  flow,  the  separation  process  is  a  function  of  the  orifice 
pressure  ratio  in  a  compressible  flow.  The  level  of  separation  decreases  as 
the  pressure  ratio,  r,  decreases .  Although  the  expanded  vena  contracta  area 
reduces  the  extent  of  the  separated  flow,  the  orifice  losses  with  compressible 
flow  are  not  decreased. 

The  most  familiar  characteri Stic  of  compressibility  associated  with  orifice 

performance  is  shown  in  Figure  24.  The  flow  through  orifices  differs  from 

that  through  most  restrictions  in  the  supercritical  behavior.  When  a  typical 

restriction  becomes  critical  (l“l  =  1.0  at  some  flow  location),  the  flow 

parameter,  <p,  is  maximized  with  respect  to  pressure  ratio.  The  restriction 

is  said  to  be  choked  at  that  location.  Then  the  adiabatic  flow  rate  becomes 

linearly  dependent  upon  upstream  total  pressure,  but  independent  of  further 

reductions  in  restriction  pressure  ratio,  r.  On  the  contrary,  however,  an 

orifice  can  be  critical  at  its  vena  contracta  (M  =  1)  and  still  exhibit  an 

c 

increasing  flow  parameter,  0,  as  the  pressure  ratio,  r,  is  reduced  This 
behavior  results  from  the  influence  that  the  downstream  static  pressure  exerts 
on  the  "free-jet"  vena  contracta.  In  highly  separated  flows  like  those 
encountered  with  sharp— edged  contractions  and  flow  angularity,  these 
compressibility  effects  become  exaggerated  for  thin  (short  duct)  orifices 
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FIGURE  24  FLOW  THROUGH  SHARP-EDGED  ORIFICES  COMPARED  TO 
FLOW  THROUGH  AN  IDEAL  NOZZLE.  (REFERENCE  (29 
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for  thicker  or  if  k  »‘i  i  ri  coaxial  flow,  the  veria  cnritraita  oriy  i  r  oriniorit  i  ■• 
isolated  from  the  downstream  static  conditions.  For  short  holes  with  an 
fi./  0  >  4  choking  will  rarely  occur  at  the  vena  contracta.  In  any  event, 
thick  orifices  and  relatively  large  orifices,  where  small  area  changes  reduce 
the  separation,  are  less  sensitive  to  the  supercritical  compressibility 
effects.  Small  edge-breaks  or  leading  edge  radii  radically  suppress 
separation.  The  geometry  of  an  orifice: 

sudden  contraction  ratio,  ft, /A  =  A 

1  o 

inlet  flow  angularity, 

leading  edge  sharpness, 

orifice  thickness,  E./HD 

sudden  expansion  ratio,  A^/A^  =  K 

affects  orifice  performance  in  either  incompressible  or  compressible  flows 

However,  the  orifice  pressure  ratio  only  influences  the  performance  of 

4 

orifices  in  compressible  flow  .  The  greatest  effect  is  seen  on  thin,  static 
orifices  operating  in  or  near  the  supercritical  regime. 

Synthesis  of  an  Orifice  Model 

Orifices  encountered  in  gas  turbine  internal  flow  systems  can  encompass  any 
combination  of  geometric  variables  important  to  flow  capacity.  Figure  25 
Therefore,  a  comprehensive  model  for  orifice  flow  characteri stic  prediction  is 
required.  The  model  proposed  by  Dodge  (24)  is  amenable  to  modifications  to 
achieve  this  flexibility: 

ko  =  "sc  +  ke  ♦  ,f  +  k5e 

Sudden  Contraction 


k 

sc 


C  C 
a  r 


# 

k 

sc 


[l-A] 


where  flow  angularity  correction  is  provided  by  C  -  ]  (0,  flow 

a 

Figure  20,  and  leading  edge  sharpness  correction  is  provided  by 
for  chamfers  or  Cr  =  z  (r/  0)  for  radii  from  Figure  21 


d l roc  t l on ) 

C  =  y  (a,  S. 
r  sc 


4Cavitating  effects  can  be  important  in  liquid  flows 
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The  sudden  contraction  model  for  orifice  applications  is  based  on  the 
experimental  work  of  Weisbach  testing  short  orifices  and  Freeman  testing 
nozzles  with  water  Both  of  these  invest igators  obtained  contraction 
coefficient  and  discharge  coefficient  data  which  resulted  in 

« 

k  =05  for  sma 11  A . 

sc 

Notice  that  the  sudden  contraction  model  is  not  the  same  model  selected  from 
reference  (26)  The  sudden  contraction  data  acquired  by  Benedict,  Carlucci, 
and  Swetz  were  for  entrances  to  long  ducts  where  separation  reattachment  and 
velocity  profile  recovery  were  attained. 

Vena  Contracta 

The  "long  hole"  k-factor  employed  by  Dodge  (24)  corrects  for  the  interactions 
of  contraction  and  expansion  geometry  with  the  incompressible  free-jet  vena 
contracta 

1/2  *  1/2 
kfi.  =  (CaCr)  k*  [1-A] 

# 

where  k ^  is  given  for  a  static  orifice  by  Figure  23, 

Sudden  [ xpans ion 


he  Sudden  expansion  mode!  of  Carnot- Borda  is  used  as  recommended  for 
ns ompre s s ! b 1 e  flow 


*  2 


;  r 


a'  ior,  t.‘  t  tie  sudden  expansion  k-factor  to  the  impact  pressure,  P-p, 
i '  i  ’ e  exit.  rather  than  the  d  y nam ic  pressure  recommended  for  the 
•'a  ', nr  ana  jpn™  contracta  losses,  generalizes  the  expansion 
*  r’*'  c  on  p»  a  •  s  :  r. .  e  f  iow  r*,qime  This  appro*  lmat  ion  is  very  good  for 


Although  the  Dodge  (24)  model  tor  orifice  c har ar ter  1 s 1 1 c s  is  specifically  for 
incompressible  flow,  it  can  be  applied  to  certain  thick  and/or  large  hole 
orifices  in  compressible  flow  with  good  accuracy.  A  technique  has  been 
developed  to  adapt  it  to  thin,  small  orifices  where  compressibility  effects 
are  most  pronounced . 


Orifice  Characteristics  for  Compressible  Flow 

The  thin  plate,  fi,/  J  <  0.1,  static  orifice,  A  =  A.  <  0.1,  is  the  most 
familiar  category  that  exhibits  strong  compressible  flow  effects  near  and  in 
the  supercritical  regime.  It  can  be  demonstrated  numerically,  however,  that 
small  blunting  of  the  leading  edge,  slight  lengthening  of  the  hole  >  : 

or  restricting  the  contraction  and/or  expansion  ratio  can  mitigate  the 
compressible  vena  contracta  characten stic  quite  rapidly.  Consequently,  a 
relatively  limited  range  of  orifice  geometry  requires  a  more  sophisticated 
compressible  flow  analysis  than  that  provided  by  the  Dodge  model.  The 
extension  of  the  Dodge  model  to  these  special  orifice  flows  and  the  range  of 
applicability  will  be  discussed. 


Perry  (29)  demonstrated  that  highly  separated  compressible  flow  through 
thin-plate,  static  orifices  behaves  linearly  in  the  subcritical  region 


<D  =Vm  ( l-r  ) 


when  modeled  in  elliptical  coordinates 


$>  =  m 


where  m  =  0.216  for  air.  Perry  (29)  represented  the  supercritical  region  as 


<t>  =  (a  +  nr)  V  l-r 


•SB 


which  in  elliptical  coordinates  becomes 


Otl 


$  =  (a  +  n/l-fD2  (1  -  V" 1-p 


where  a  =  0.449  and  n  =  0.241  for  air. 


The  slope  of  the  supercrit ical  flow  model  was  found  to  be  88%  of  the 
subcritical  model  slope  at  the  choking  pressure  ratio, 

* 

m  =  0.88  m 
c 

The  compressible  flow  model  proposed  by  Perry  (29)  for  thin,  static  orifices 
can  be  generalized  to  model  the  limited  range  of  orifices  where  vena  contracta 
compressibility  is  important.  Table  VI.  The  orifice  model  by  Dodge  can  be 
used  to  generate  the  incompressible  flow  characten sties  for  the  specific 
orifice  geometry  des,red.~>  A  common  slope  for  the  compressible  and  the 
incompressible  flow  can  be  found  near  r  =  0.87  as 


m  =  (I') 

>r  ~  0.87 

n  =\GT  /-L-fi-  0.88/j±-\i\ 

I  l-r*  ll-r*  [_  2  ^  1+r*  /  J  J 


a  “  nr* 


where 


=  (vTl) 


Y/(Y-1) 


for  Y  =  1 . 4 


*  2 
$  =  m  [l-(r*)  ] 


&A  flow  characteristic  curve  generator  similar  to  DU1  in  BC88  PLUS  can  be 
used  to  solve  for  $  at  r  using  the  k-factors  from  the  Dodge  orifice  model 


»Vi  yt 
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Table  VI . 

Applicable  range  of  the  compressible  flow  parameters 
for  orifice  models  based  on  Perry  (29). 

k  s  c  *  0.4 

k£  >  0.65 

A.  <  0.125  which  is  the  same  as  (k+  >  0.765) 

m  <  0.26  SG 

It  can  be  seen  from  the  supercritical  orifice  model  that  a  =  <t>  when  r  =  0. 
Therefore,  a  >  0  532  for  air  cannot  be  allowed. 

Other  more  stringent  flow  modeling  restrictions  place  tighter  limits  on  the 
applicable  range  of  parameters  The  entrance  flow  must  be  severely  separted, 
and  the  orifice  hole  must  be  short  enough  to  preclude  flow  reattachment  within 
the  hole.  The  orifice  exit  area  ratio  must  provide  a  large  expansion  so  that 
the  flow  reattaches  in  the  far  downstream  field  of  the  tangent  duct.  If  these 
stipulations  as  quantified  in  Table  V  are  met,  the  compressible  flow 
characteristics  of  the  orifice  can  generally  be  modeled  satisfactorily  over 
the  complete  range  of  pressure  ratio  as  outlined. 

Example  Calculations  for  Generalized  Orifice  Flow  Characteristics 

Two  orifices  have  een  modeled  to  demonstrate  the  calculation  procedures  for 

1)  A  conventional  thick-plate  orifice  with  nozzle-like  characten sti cs 

2)  A  generalized  thin  orifice  with  definite  vena  contracta  compressibility 
characteri st i c  s 

The  detailed  calculations  for  restriction  1  and  restriction  2  are  located  in 
the  Derivations  section  of  the  Appendix  The  Allison  Gas  Turbine  Engines 
version  of  a  flow  characteri Stic  curve  generator  program  titled  DU1  was 
utilized  to  calculate  the  airflow  parameter,  5>,  as  a  function  of  the  total 
pressure  ratio,  P^/P^  Constant  values  for  restriction  k-factors  were 
used  with  the  exception  of  the  E ANNO  wall  friction  calculation  The  Moody 
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correlation  for  Fanning  friction  factor  was  calculated  at  each  flow  condition 
on  the  curve 


Restriction  1  is  nozzle-like  according  to  the  low  internal  k-factors  The 
flow  charac ter  1 s t i c  curve  generated  by  DU1  calculations  of  the  Dodge  model 
represents  the  thick-plate  orifice  performance  very  well.  This  orifice  will 
exhibit  a  classical  choked  flow  charac ter  1 s t i c  because  of  the  minimal  entrance 
separation  and  the  internal  flow  reattachment. 


Restriction  2  exhibits  definite  thin  orifice-like  performance  according  to  the 
combination  of  high  k-factors  for  the  loss  elements  and  relatively  low  slope' 
in  the  elliptical  parameters  Since  sustained  entrance  flow  separation  with 
free  jet  vena  contracta  characteri st ics  are  indicated,  the  flow  curve  will  not 
choke  but  will  continue  to  rise  in  the  superc r l t l ca 1  regime  Therefore,  the 
Perry  model  was  used  to  predict  the  orifice  flow  curve.  The  Dodge  model  was 
calculated  in  the  DU1  program  to  determine  the  orifice  baseline  performance  in 
the  low  pressure  ratio,  P^/P^,  or  "incompressible"  regime. 


The  orifice  flow  characteristics  are  derived  directly  from  the  component  loss 
k-factors  in  the  DU1  calculation.  However,  if  a  curve  of  k-factor  versus  flow 
parameter  is  required  for  inserting  an  orifice  into  a  component  restriction  of 
a  more  extensive  geometry,  the  "kurve"  data  can  be  developed  from  the  flow 
curve  at  each  P^/P^  as  follows 


<t 


k  _ 


1  -  (P  /P  ) 
U  U 


(see  note) 


The  "ku’-ve  1"  l  n  tn<>  example  DUi  input  is  the  Perry  static  orifice  curve  o* 
x  factor  versus  flow  ‘actor.  $  7 h i s  data  file  contains  all  of  the  element 

"k  turn 


for  f  i  r  i  f  i  •  e  modeling  usi  ng  t  he  pc  oc  edure  s  o*  [lodge  ,  Per  r  y  or 


The  calculation  of  (p)  from  ®  requires  the  solution  of  the  implicit  equati 


..V 


>g 


\  -J±1- 
1  *  2=1  M2  2(V-1) 


for  Mach  number.  The  (  p  )  can  be  found  directly  from 


(?)  -7 


M 


p  /  2  y 

(l+^i„2)>TT 
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TOTAL  NUMBER  OF  REFERENCES  FOR  LITERATURE  SURVEYS  =  4 


NOMENCLATURE 


et  Small  or  minor  dimension  of  a  cross-section  of  a  rectangular  or 

elliptical  duct,  in 

a  Limiting  9  of  an  orifice  at  r  =  0 

^  Area  ratio,  An/A3 

A  Cross-sectional  area  of  a  duct,  in  .2 

Are^  Real  or  defined  reference  area  for  a  flow  restriction,  in  ? 

b  Large  or  major  dimension  of  a  cross-secton  of  a  rectangular  or 

elliptical  duct ,  in . 

Cp  Specific  heat  of  the  air  at  constant  pressure,  Btu/lbm  deg  R 

cu  Specific  heat  of  the  air  at  constant  volume,  Btu/lbm  deg  R 

Ca  Angularity  correction  for  obliquely  incident  flow  into  a  restriction 

Cft  Area  change  correction  for  bends 

Cc  Contraction  coefficient,  Ac/An 

C(j  Drag  coefficient 

Cq  Discharge  coefficient  for  a  restriction,  m/micj 

Cf  Frictional  influence  coefficient  for  aerodynamical ly  rough  wall 

Cg  Cros s-sect lona 1  area  geometry  correction 

Cj  Influence  coefficient  for  a  general  parameter  which  is  different  than 

that  for  the  reference  restriction 

C£  Downstream  tangent  correction 

Compressibility  influence  coefficient  for  high  velocity  flow 

Cr  Edge  break  correction  for  a  restriction  area  reduction 

Cu  Velocity  coefficient  for  a  restriction,  V/Vjd 

d  Diameter  of  a  circular  cross-section,  or  small  diameter  of  an  annular 

cross-sectional  area  of  a  duct,  in. 

D  Large  diameter  of  an  annular  cross-sectional  area  of  a  duct,  in 

f  Fanning  friction  factor  for  flow  in  straight  ducts 
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fc 

f(  ) 
fd 

9c 

h 


D 

k 

k+ 


NOMENCLATURE  (con't) 

Darcy-Wei sback  friction  factor  for  flow  in  curved  ducts 
Functional  relationship  of  independent  variables  (  ) 

Aerodynamic  drag  force,  lbf 

Conversion  factor,  32.174  lbm  ft/lbf  sec 

Height  (maximum)  of  a  cross-section  of  a  duct  in  the  plane  defined  by 
radius  r  normal  to  the  bend  axis  (h  can  be  the  same  as  a  or  b  of  a 
rectangular  or  elliptical  duct  for  example),  in 

Hydraulic  diameter  of  a  cross-section  area  of  a  duct,  in. 

Total  pressure  loss  coefficient  based  on  q 

Total  pressure  loss  coefficient  based  on  (P  -  p) 

Total  pressure  loss  coefficient  for  the  reference  restriction 
(usually  in  the  incompressible  flow  regime) 


fi.  Length  along  the  centerline  of  a  duct,  in. 
in  Subscript  slope  of  an  orifice  performance 
fti  Mass  flowrate,  lbm/sec 


ifild  Ideal  m  which  would  pass  through  a  lossless  restriction  if  the 

avialable  cross-sectional  area  flowed  full,  lbm/sec 


M 

n 

P 

P 

(P  -  P) 

r 

q 

« 

9 


Mach  number 

Supercritical  constant  for  orifice  performance 

Static  pressure  of  the  air,  psia 

Total  pressure  of  the  air,  psia 

Impact  pressure  of  the  air,  psia 

Perimeter  of  a  duct  cross-sectional  area,  in. 

Dynamic  pressure  of  the  air,  psi 

Volumetric  or  mass  flow  ratio,  Qn/Q3  or  mn/m3 

Volume  flowrate,  ft^/sec 

Radius  of  curvature  for  the  centerline  of  a  ciruclar-arc  bene,  or 
edge  break  or  fillet  radii  at  tube-wall  intersections,  in. 


* 

< 

< 

4 

i 

3 
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NOMENCLATURE  (con't) 


Orifice  pressure  ratio,  PQ/PU 

Radius  of  a  circular  cross-section,  in. 

Relative  radius,  r/  D 

Elliptical  pressure  ratio  function  for  orifices,  l-r^ 

Specific  gas  constant  of  the  air,  lbf  ft/lbm  deg  R 
Reynolds  number 

Total  temperature  of  the  air,  deg  R 
Velocity  of  the  air,  ft/sec 

Velocity  of  an  isentropie  one-dimensional  flow  filling  the  same  area, 
ft/ sec 

Cartesian  coordinate  or  an  arbitrary  geometrical  variable 
Cartesian  coordinate  or  an  arbitrary  geometrical  variable 
Turning  loss  term  for  bend  k-factor  equations 
Complementary  bend  angle,  degrees 

Ratio  of  specific  heats  of  the  air,  cp/cu 
Effective  "sand  grain"  wall  roughness,  in. 

Bend  angle,  degrees 

Sudden  expansion  area  ratio,  A1/A2 

Sudden  contraction  area  ratio,  A2/A1,  or  Aj/Aq  fpr  profoces 
Dynamic  viscosity  of  the  air,  lbm/ft  sec 
Static  density  of  the  air,  lbm/ft^ 

Compressible  flow  parameter,  m  -/T/P  A,  lbm  °R1^2/lbf  sec 
Elliptical  flow  parameter  function  for  orifices,  <t>2 
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NOMENCLATURE  (con't) 


Superscripts 


Average  or  effective  value  of  a  nonconstant  parameter 


Critical  condition,  M  =  1.0 


Subscripts 


Bend  restriction 

Fluid  stream  contraction  due  to  separation  at  an  abrupt  flow  area 
reduction,  primarily  a  vena  contracta 

Downstream  tangent  duct 

Downstream  of  a  restriction 

Restriction  exit  area  component  (discharge  contribution) 

Inside  wall 

General  location  in  the  internal  flow  system 
Outside  wall 

Sudden  contraction  restriction 
Sudden  expansion  restriction 
Upstream  tangent  duct 
Upstream  of  a  restriction 
Velocity  of  the  fluid  stream 

Free  stream  condition  upstream  of  a  flow  obstacle 
Inlet  aea  of  a  flow  restriction 
Exit  area  of  a  flow  restriction 


Junction  or  branch  leg  carrying  the  combined  flow 
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APPENDIX 


Summary  of  Derivations 

This  Appendix  contains  the  detailed  calculations  for: 

Restriction  l--0r1f1ce  flow  characteristics  for  a  nozzle-like  geometry.  The 
Dodge  (24)  model  was  used  to  obtain  k-factors  for  the  DU1  program. 

Restriction  2--0r1f1ce  flow  characteristics  for  a  thln-plate  geometry.  The 
Perry  (29)  model  was  used  to  develop  the  flow  characteristic  curve.  The  Dodge 
(24)  model  was  used  to  obtain  k-factors  for  the  DU1  program. 

The  flow  characteristic  curve  was  converted  to  an  overall  k-factor  curve 
(KURVE2)  suitable  for  Internal  Dill  use,  as  demonstrated  by  the  Restriction  7 
calculation. 

Restriction  3 — The  Dodge  model  was  used  to  generate  an  orifice  flow 
characteristic  representative  of  the  minimum  component  losses  of  Table  VI.  An 
"Incompressible"  flow  slope  of  m  =  0.317  was  found  at  r  ~  0.87.  This  yields 
a  value  of  a  =  0.544  which  exceeds  the  theoretical  limit  of  a  =  0.532  for 
air.  Therefore,  the  need  to  observe  the  overall  slope  limit,  m  <  0.26,  In 
addition  to  the  component  loss  limits  Is  demonstrated. 

Restriction  4--Prov1des  a  comparison  of  the  static  orifice  curve  calculated  by 
the  Dodge  model  with  the  test  data  correlation  of  Perry  (29).  The  need  for  an 
alternative  model  for  highly  separated  orifices  can  be  seen  from  the  flow 
characteristic  discrepance  at  high  Py/Pg- 

Restriction  5--Demonstrates  the  use  of  a  k-factor  curve  for  the  accurate 
representation  of  orifice  performance  In  a  DU1  restriction  calculation.  The 
CURVE1  example  Is  for  the  static  orifice  of  Perry  (29). 

Restriction  6--Demonstrates  the  rapid  expansion  of  the  vena  contracta  In 
high-speed  compressible  flow.  It  can  be  seen  that  1/  D  *  1.2  Is  the 
theoretical  limit  for  a  Perry  model  of  a  static  orifice. 


Restriction  7--Ver1f1es  the  equivalence  of  the  k-factor  curve  (KURVE2)  to  the 
OU 1  flow  characteristic  calculated  for  Restriction  2. 

Restriction  8  ( 2k)--Demonstrates  an  Imbedded  k-factor  curve  (KURVE2)  for  a 

standard  component  within  a  more  estenslve  flow  characteristic  model,  W>r 

T  /P  versus  P./P.,  (A  =  constant,  as  modeled). 

U  U  U  0  n 

Restriction  8  (A*)--Demonstrates  the  generalization  of  an  extensive  flow 
characteristic  model  (In  this  case,  restriction  8(2K))  to  all  geometrically 
similar  restrictions  on  the  arbitrary  basis  of  the  minimum  flow  area,  A 

o 

( KURVE2 ) . 


4f  Moody  correlation  D  =  0.35  in. 

SL  =  0.42  in. 


A  =  A] /A0  =  0.167 
K  =  Aj/A2  =  0.033 
Sudden  contraction 


£  =  30  p  in. 

Pi  =  0  408 
P2  =  0. 183 


^sc  _  ^r  k sc C 1  —  A] 

=  (0.06)  0.5  [0.833] 
=  0.025 
Vena  contracts 


1/2  *  1/2 

k  =  [1  -  A] 1/Z  ]1  -  \] 


=  (0.245)  0.43  ]0 . 913]  [0.967] 
=  0.093 


Wall  friction 
kf  =  4f  (£/  D) 
Sudden  expansion 
k, 


>se  =  [1  ~  ^J2 
=  0.934 


(<<  0.4) 


(<<  0.65) 


Internal  calculation 


(A.  <  0.125) 


Slope 

m  =  0/r 

near  Pu/Pq  <_  1.15 

pu/po 

1  .  1283607 

1 . 1761978 

0.3351456 

0.3757693 

§ 

0.11232 

0. 14120 

0.21458 

0.27195 

m 

0.52346 

0.51922 

n 

0.37128 

0.36977 

a 

0.69829 

0 . 69545 

r* 

0.52828 

0.52828 

<J># 

0.61431 

0.61181 

A-4 


/PD 

1  1177461 

1  1710073 

4> 

0.2323073 

0.2689874 

t 

0.05397 

0.07235 

0  19959 

0  27074 

m 

0.27039  slightly  exceeds 

0  26724 

n 

0.26684 

0.26528 

a 

0.50187 

0.49894 

r* 

0,52828 

0.52828 

0.44151 

0  43893 

A-5 


Dodge  (24)  incompressible  flow  model 


:  ksc  +  k£ 

+  kf 

+  k 

se 

* 

k  =0.5 

sc 

Cr  =  0.85 

Fig. 

21 

asc/  D  =  0.02,  a  =  90 

Ca  =  1.35 

Fig. 

20 

0  =  30°  Weisbach  Eq . 

ka=  i  00 

Fig. 

23 

a/  D  =  0.5 

4f  Moody  correlation  D  =  0.35  in. 

a  0. 18  in. 

£  =  30  \i  in. 


A  =  Ai/A0  =  0.167 
A.  =  =  0.  100 

Sudden  contraction 
# 


k  =  C  C  k  ]1  -  A] 
sc  a  r  sc 

=  1.35  (0.85)  0.5  ]0 . 833 ] 

=  0.478 


01  =  0.408 
02  =  0.316 


(>  0.4) 


Vena  contracta 


ka  =  (Ca  Cr)1/2  kscCl  '  A]1/2  [1  -  W 
=  (1.071)  1.00  [0.913]  [0.900] 

=  0.880  (<  0.65) 

Wall  friction 

kf  =  4f  ( a /  D)  Internal  calculation 

Sudden  expansion 

k3e  =  [1  -  A.]2  (A.  <  0.125) 

Slope 

m  =  $/  near  Pu/Pq  <_  1 . 15 


RES 2 — SIMULATION  OF  A  GENERALIZED  ORIFICE  COMPRESSIBLE  FLOW  -  L/D  =  OS, 
EDGE-BREAK  PERRY  MOOEL  BASED  ON  DODGE  "INCOMPRESSIBLE"  PERFORMANCE  TOR  AIR 


p  /Pn 
u  D 

m  =  0.2674 
r 

8  =  0- 

2654  a  = 

Me 

o 

vO 

n  «o 

1,00 

1 . 0000 

0 

0 

0 

1.00363 

0.99638 

0.4394 

0.4788 

0.00160 

1.01780 

0.98251 

0.09629 

0. 10549 

0.00773 

1 .04161 

0.96005 

0. 14470 

0. 15989 

0.01758 

1.07493 

0.93029 

0.18968 

0.21203 

0.03050 

1.11775 

0.89465 

0.23102 

0.26189 

0.04577 

1.23561 

0.80932 

0.30373 

0.35636 

0.08142 

1 .31332 

0.76143 

0.33521 

0.40115 

0.10082 

1  . 40726 

0.71060 

0.36383 

0.44482 

0.12091 

1.52107 

0.65743 

0.38965 

0.48738 

0. 14134 

1.66159 

0.60183 

0.41297 

0. 52924 

0.16201 

1  .  86193 

0.53708 

0.43620 

0.57534 

0.18514 

1 . 89293 

0.52828 

0.43906 

0.58140 

0. 18820 

2.11115 

0.47368 

0  45326 

0.61304 

0.20413 

2.3 

0.43478 

0 . 46195 

0.63391 

0.21461 

2 . 5 

0 . 40000 

0.46881 

0.65138 

0.22334 

3.0 

0.33333 

0  47972 

0.68140 

0.23819 

3.5 

0.28571 

0.48588 

0.69985 

0.24720 

4.0 

0.25000 

0  48967 

0.71187 

0.25302 

5.0 

0 . 20000 

0  49386 

0.72585 

0.25972 

7.0 

0.14286 

0.49716 

0.73746 

0.26523 

10.0 

0. 10000 

0.49865 

0.74289 

0.26780 

20.0 

0.05000 

0.49938 

0.74560 

0.26907 

100.0 

0.01000 

0.49922 

0.74501 

0.26879 

1000.0 

0.00010 

0.49908 

0.74449 

0.26854 

k0 


2  .  255 
2 . 25726 
2 . 26281 
2 . 27240 
2 . 28549 
2.30139 
2 . 34192 
2  .  34192 
2 . 39348 
2.42366 
2.45762 
2.50037 
2 . 50654 
2.57833 
2.62267 
2 . 68652 
2 . 79890 
2 . 88946 
2  96419 
3 .08019 
3.23165 
3 . 36077 
3.35070 
3.68319 
3.72341 


//// 


////////  INPUT  DATASET  /////////////  DBA  ////////// 


KURVE 

1  20 

0.  0 

0. 0328 

0. 0b53 

0. 0925 

0. 1302 

0. 158b 

0. 1823 

0.2027 

0. 24b0 

0.2790 

0.3080 

0.3320 

0.3540 

0. 3730 

0.3890 

0.4030 

0.4230 

0. 43b0 

0. 4440 

0. 4490 

2.7990000  2. 7940058  2.7940899 

2.8042504 

2. 8195bb5 

2.8322270 

2. 84b0312 

2. 8bl 13bC  2.87b01b8  2.9434909 

2.9794882 

3. 0389848 

3. 077  b5b3 

3. 1275b57 

3. 1960778  3.2702482  3.4953703 

3. 7851911 

4. 1324977 

4. 7b88857 

KURVE 

2  20 

0.  0 

0. 04394 

0. 09b29 

0. 14470 

0. 189b8 

0.23102 

0. 2b907 

0.30373 

0.33521 

0. 3b383 

0. 389b5 

0. 41297 

0.4390b 

0.4532b 

0. 4bl95 

0. 4b881 

0.47972 

L. 48588 

0. 4938b 

0.49908 

2. 25500 

2.2572b 

2. 2b28 1 

2.27240 

2.28549 

2.30139 

2. 32019 

2. 34192 

2. 3bb23 

2. 39348 

2. 423bb 

2. 4S7b2 

2. 50b54 

2.57833 

2. b33b7 

2. b8bS2 

9 

2. 79890 

2.8894b 

3. 08019 

3. 72341 

RES 

1 

SIMULATION  OF  A  GENERALIZED  ORIFICE  COMPRESSIBLE  FLOW  24APR8b 

BASED  ON 

THE  DODGE  MODEL  -  L'D 

=  1.2,  LEADING  EDGE 

RADIUS 

5 

0.  10 

1.40  28.97 

PLOT 

0.  bO 

540. 

0.  10 

540. 

0. 02S 

Q 

0.  10 

540. 

0.  093 

Q 

0.  10 

540. 

0.42 

0.35  30. 

3.  00 

540. 

125. 

2 

0.934 

PT 

RES 

2 

SIMULATION  OF  A  GENERALIZED  ORIFICE  COMPRESSIBLE  FLOW  24APR8b 

BASED  ON 

THE  DODGE  MODEL  -  L/D 

=  .5,  EDGE  BREAK,  OBLIQUE  ANGLE 

5 

0. 10 

1.40  28.97 

PLOT 

0.  bO 

540. 

0.  10 

540. 

0.478 

Q 

0.  10 

540. 

0. 880 

Q 

0. 10 

540. 

0.  18 

0.35  30. 

1.  00 

540. 

125. 

3 

0.810 

PT 

RES 

3 

SIMULATION  OF  A  COMPRESSIBLE  FLOW  STATIC 

ORIFICE 

GFH  24APRBb 

LIMIT 

OF  THE  DODGE  MODEL  — 

TABLE  VI  COMPONENT 

LOSSES 

5 

1.  0 

1.40  28.97 

PLOT 

1000. 

100. 

1.  00 

100. 

0.40 

Q 

1.00 

100. 

0.  b5 

Q 

1.  00 

100. 

0.  0 

1000. 

100. 

10.0 

0 . 7bb 

PT 

4 


RES  6  SIMULATION  OF  A  COMPRESSIBLE  FLOW  STATIC  ORIFICE  GFH  24APR86 
BASED  ON  PERRY  EMPIRICAL  DATA  AND  DODGE  INCOMPRESSIBLE  MODEL 


5 

1.  0 

1.40  28.97 

PLOT 

1000. 

100. 

1.  00 

100. 

0.50 

Q 

1.  00 

100. 

1.34 

Q 

1.00 

100. 

0.  0 

1000. 

100. 

10.0 

5 

1.0 

PT 

RES  5 

STATIC  ORIFICE  RESTRICTION 

IN  COMPRESSIBLE  FLOW 

GFH 

24APR8b 

BASED  ON 

THE  K-FACTOR  DATA 

TABLE 

FOR  EMPIRICAL  PERRY 

MODEL 

4 

1.  0 

1.40  28.97 

PLOT 

10.000 

100. 

1. 0000 

100. 

1.0000 

100. 

1  Q 

10. 000 

100. 

10. 

b 

RES  fa 

SIMULATION  OF  A  COMPRESSIBLE  FLOW 

STATIC  ORIFICE 

GFH 

24APR8b 

LIMIT  OF 

THE  DODGE  MODEL  - 

L/D 

=  1.2 

5 

1.  0 

1.40  28.97 

PLOT 

1000. 

100. 

1.00 

100. 

0.50 

Q 

1.  00 

100. 

0.  43 

Q 

1.  00 

100. 

0.024 

A-9 


1000. 

100. 
10. 0 

7 

1.  0 

PT 

RES 

7( 2K  ) SIMULATION  OF  A  GENERALIZED  ORIFICE 

COMPRESSIBLE  FLOW  24APR8b 

BASED 

ON 

THE  DODGE  MODEL  -  L/D  =  .5. 

EDGE 

BREAK.  0BLI9UE  ANGLE 

4 

0. 10 

1.40  28.97 

PLOT 

0.  bO 

540. 

0.  10 

540. 

0. 10 

540. 

2 

Q 

1.  00 

540. 

125. 

8 

RES 

8 ( 2K  JCOMPLEX 

RESTRICTION  MODEL  INCLUDING 

GENERALIZED  ORIFICE  RES  2 

BASED 

ON 

THE  DODGE  MODEL  -  L/D  =  .5, 

EDGE 

BREAK.  OBLIQUE  ANGLE 

8 

1.  00 

1.40  28.97 

PLOT 

1000. 

bOO. 

1. 00 

bOO. 

0. 578 

Q 

1. 00 

580. 

5.50  1.12850. 

0.  bO 

540. 

0. 050 

Q 

0. 10 

540. 

2 

Q 

1. 00 

540. 

1.  00 

520. 

3.00  1.12850. 

1000. 

500. 

125. 

9 

1.  0 

PT 

RES 

8( A» ICOMPLEX 

RESTRICTION  MODEL  INCLUDING 

GENERALIZED  ORIFICE  RES  2 

BASED 

ON 

THE  DODGE  MODEL  -  L/D  =  .5. 

EDGE 

BREAK.  OBLIQUE  ANGLE 

8 

0.  10 

1.40  28.97 

PLOT 

1000. 

bOO. 

1. 00 

bOO. 

0.578 

9 

1.  00 

580. 

5.50  1.12850. 

0.  bO 

540. 

0.  050 

Q 

0. 10 

540. 

2 

9 

1.00 

540. 

1.  00 

520. 

3.00  1.12850. 

1000. 

500. 

125. 

1.  0 

PT 
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